Science China Life Sciences

Leuciscus waleckii is widely distributed in Northeast Asia and has high economic value. Different from its freshwater counterparts, the population in Lake Dali Nur has a strong alkalinity tolerance and can adapt to extremely alkaline-saline water with bicarbonate over 50 mmol/L (pH 9.6), thus providing an exceptional model with which to explore the mechanisms of adaptive evolution under extreme alkaline environments. Here, we assembled a high quilty chromosome-level reference genome for L. waleckii from Lake Dali Nur, which provides an important genomic resource for the exploitation of alkaline water fishery resources and adaptive evolution research across teleost fish. Notably, we identified significantly expanded long terminal repeats (LTRs) and long interspersed nuclear elements (LINEs) in L. waleckii compared to other Cypriniformes fish, suggesting their more recent insertion into the L. waleckii genome. We also identified expansions in genes encoding gamma-glutamyltransferase, which possibly underlie the adaptation to extreme environmental stress. Based on the resequencing of 85 L.waleckii individuals from divergent populations, the historical population size of L.waleckii in Lake Dali Nur dramatically expanded in a thousand years approximately 13,000 years ago, and experienced a cliff recession in the process of adapting to the alkaline environment of Lake Dali Nur approximately 6,000 years ago. Genome scans further revealed the significant selective sweep regions from Lake Dali Nur, which harbour a set of candidate genes involved in hypoxia tolerance, ion transport, acid-base regulation and nitrogen metabolism. In particular, 5 alkali population specific nonsynonymous mutations were identified in CA15 gene copies. In addition, two sites with convergent amino acid mutation were detected in the RHCG-a gene among several alkali environment adapted Cypriniformes fish, this mutation may increase the NH3 excretion rate of the RHCG channel. Our findings provide comprehensive insight into the genomic mechanisms of L. waleckii and reveal their adaptative evolution under extreme alkaline environments.

Cloned animal have been reported by somatic cell nuclear transfer (SCNT) for many years. However, the SCNT is extremely inefficient, and zygotic genome activation (ZGA) is required for SCNT and chemical mediated reprogramming. To identify candidate factors that facilitate ZGA in reprogramming, we performed siRNA-repressor and mRNA-inducer screening, which revealed Dux, Dppa2, and Dppa4 as key factors enhancing the ZGA in SCNT. Direct injection of ZGA-inducers had no significant effect on the SCNT blastocyst formation, and even destroyed the ZGA. Through a inducible Dux transgenic mouse model, we demonstrate that transient overexpression of Dux not only improved SCNT efficiency, but also increased the efficiency of chemical reprogramming. Transcriptome profiling revealed that Dux treated SCNT embryos were similar to fertilized embryos. Furthermore, transient overexpression of Dux combined with inactivation of DNA methyltransferases (Dnmts) further promote the overall development of SCNT-derived animals. These findings enhance our understanding of ZGA-regulators in somatic reprogramming.

Functional profiling on whole-metagenome shotgun sequencing (WMS) has made great contribution to the development of our understanding in microbe-host interactions. In this work, we revealed that severe microbial functional information loss of current functional profiling methods existed at both taxon-level and community-level. To correct the distortion brought by information incompleteness, we developed a new framework, RFW (Reference based functional profile inference on WMS), to infer microbial functional abundance on WMS through utilizing information from genome function annotation and WMS taxonomic profile. Furthermore, we built up a new algorithm for absolute abundance change quantification of microbial function between groups under RFW framework. By applying RFW to several datasets related to autism spectrum disorder and colorectal cancer, we revealed that RFW greatly renewed our knowledge in downstream analysis, including differential microbial function identification, association analysis between microbial function and host phenotype, etc. RFW are open-source and freely available at https://github.com/Xingyinliu-Lab/RFW.

As members of the -proteobacteria group, Caulobacter crescentus and its relatives are known for their asymmetric life cycle and comprehensive applications in gene delivery, agricultural biotechnology, and the production of high-value compounds. However, genetic manipulations of these bacteria are often time-consuming and labor-intensive due to the lack of efficient genome editing tools. Here, we report a practical CRISPR/SpCas9M-reporting system that overcomes the limitations of SpCas9 expression, enabling efficient, markerless, and rapid genome editing in C. crescentus. As a demonstration, we successfully knocked out two genes encoding the scaffold proteins, achieving apparent editing efficiencies up to 80%. Key components, including the Cas protein, Cas inducer, sgRNA, homologous arms, and reporter, were systematically analyzed and optimized to enhance the editing efficiency or decrease the cell lethality. A nearly zero off-target ratio was observed after the curing of the editor plasmid in editing strains. Furthermore, we applied the CRISPR/SpCas9M-reporting system to two C. crescentus relatives, Agrobacterium fabrum and Sinorhizobium meliloti, establishing it as an efficient and reliable editing strategy. We anticipate that this system could be applied to other hard-to-edit organisms, accelerating both basic and applied research in -proteobacteria.

The COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has overwhelmed many health systems globally. Here, we aim to identify biological markers and associated biological processes of COVID-19 using a bioinformatics approach to elucidate their potential pathogenesis. The gene expression profile of the GSE152418 dataset was originally produced by using the high-throughput Illumina NovaSeq 6000. Kyoto Encyclopedia of Genes and Genomes pathway (KEGG) and Gene Ontology (GO) enrichment analyses were applied to identify functional categories and biochemical pathways. KEGG and GO results suggested that biological pathways such as "Cancer pathways" and "Insulin pathways" were mostly affected in the development of COVID-19. Moreover, we identified several genes including EP300, CREBBP, and POLR2A were involved in the virus activities in COVID-19 patients. We further predicted that some inhibitors may have the potential to block the SARS-CoV-2 infection based on the L1000FWD analysis. Therefore, our study provides further insights into the underlying pathogenesis of COVID-19.

We tracked the effective reproduction number Rt of SARS-CoV-2 Omicron BF.7 in Beijing in November - December 2022 by fitting a transmission dynamic model parameterized with real-time mobility data to (i) the daily number of new symptomatic cases on November 1-11 (when the zero-covid interventions were still strictly enforced) and (ii) the proportion of individuals who participated in online polls on December 10-22 and self-reported to have been previously test-positive since November 1. After the announcement of "20 measures", we estimated that Rt increased to 3.44 (95% CrI: 2.82 - 4.14) on November 18 and the infection incidence peaked on December 11. The cumulative infection attack rate (i.e. the proportion of population who have been infected since November 1) was 43.1% (95% CrI: 25.6 - 60.9) on December 14 and 75.7% (95% CrI: 60.7 - 84.4) on December 22. Surveillance programmes should be rapidly set up to monitor the evolving epidemiology and evolution of SARS-CoV-2 across China.

As a promising biotechnology, fish germ cell transplantation shows potentials in conservation germplasm resource, propagation of elite species, and generation of transgenic individuals. In this study, we successfully transplanted the Japanese flounder (P. olivaceus), summer flounder (P. dentatus), and turbot (S. maximus) spermatogonia into triploid Japanese flounder larvae, and achieved high transplantation efficiency of 100%, 75-95% and 33-50% by fluorescence tracking and molecular analysis, respectively. Eventually, donor-derived spermatozoa produced offspring by artificial insemination. We only found male and intersex chimeras in inter-family transplantations, while male and female chimeras in both intra-species and intra-genus transplantations. Moreover, the intersex chimeras could mature and produce turbot functional spermatozoa. We firstly realized inter-family transplantation in marine fish species. These results demonstrated successful spermatogonial stem cells transplantation within Pleuronectiformes, suggesting the germ cells migration, incorporation and maturation within order were conserved across a wide range of teleost species.

As reported previously, exosomes have significant impacts on the physiological and pathological state in vivo. Exosomes have been extensively studied as a drug delivery carrier and some exosomes drugs have been undergoing clinical research. The mechanisms of exosome production, transport, and secretion remain to be studied in depth. Thus, here, we developed a novel CRISPR-UMI-based gene screening system, with which we can reveal the exosomes biogenesis, transport, and uptake mechanisms on a genome-wide scale. This system consists of two parts: one is a gene knockout component; the other part is a unique molecular identifiers (UMI) labeling component that can label the exosomes produced by the knockout cells, in which each sgRNA corresponds to a specific UMI. In this way, by detecting the UMI loaded in the exosomes, we can trace the knockout gene. In this study, we first verified the function of each component using plasmids and lentiviruses respectively. Next, we simulated the infection of cells with lentiviral libraries using a single lentivirus to validate the functionality of the screening system. Finally, we constructed a CRISPR-UMI-based library targeting 15 genes (genes with clear effects on exosomes biogenesis) to further validate the reliability of the screening system. The development of this screening system is of indispensable importance for the in-depth study of the mechanisms of exosome production and secretion, as well as for the improvement of exosomes production and the advancement of exosomes industrialization.

Corals create an ecosystem, called a holobiont, with intracellular algae (zooxanthellae) and resident bacteria. Zooxanthellae and some bacteria play major roles in the physiological properties of the coral host. However, because of the difficulties in experimental verification of cross-organism interactions, the mechanisms underpinning these interactions are largely unknown. To address this, we here generated and then analyzed multi-omics datasets for corals, zooxanthellae, and bacteria collected at Okinawa, Japan, from November 2014 to September 2016. Using cross-organism co-expression analysis, we successfully characterized the host-alga relationship in the coral holobiont. Specifically, we observed that the coral host dominates the zooxanthellae. The multi-omics analysis also suggested that infection with coral-associated bacteria Endozoicomonas likely involves coral-like ephrin ligands, triggering an immune response of the coral host. This study highlights the potential of the multi-omics approach to elucidate coral-microbe interactions.

Candidatus Parvarchaeales, representing a DPANN archaeal group with limited metabolic potentials and reliance on hosts for their growth, were initially found in acid mine drainage (AMD). Due to the lack of representatives, however, their ecological roles and adaptation to extreme habitats such as AMD, as well as how they diverge across the lineage remain largely unexplored. By applying genome-resolved metagenomics, 28 Parvarchaeales-associated metagenome-assembled genomes (MAGs) representing two orders and five genera were recovered. Among them, we identified three new genera and proposed the names Candidatus Jingweiarchaeum, Candidatus Haiyanarchaeum, and Candidatus Rehaiarchaeum with the former two belonging to a new order Candidatus Jingweiarchaeales. Further analyses of metabolic potentials revealed substantial niche differentiation between Jingweiarchaeales and Parvarchaeales. Jingweiarchaeales may rely on fermentation, salvage pathways, partial glycolysis, and pentose phosphate pathway (PPP) for energy reservation, while the metabolic potentials of Parvarchaeales might be more versatile. Comparative genomic analyses suggested that Jingweiarchaeales are more favorable to habitats with higher temperatures and Parvarchaeales are better adapted to acidic environments. We further revealed that the thermal adaptation of these lineages especially for Haiyanarchaeum might rely on innate genomic features such as the usage of specific amino acids, genome streamlining, and hyperthermal featured genes such as rgy. Notably, the acidic adaptation of Parvarchaeales was possibly driven by horizontal gene transfer (HGT). Reconstruction of ancestral states demonstrated that both may originate from thermal and neutral environments and later spread to mesothermal and acidic environments. These evolutionary processes may also be accompanied by adaptation toward oxygen-rich environments via HGT. ImportanceCandidatus Parvarchaeales may represent a lineage uniquely distributed in extreme environments such as AMD and hot springs. However, little is known about the strategies and processes of how they adapted to these extreme environments. By the discovery of potential new order-level lineages - Jingweiarchaeales and in-depth comparative genomic analysis, we unveiled the functional differentiation of these lineages. Further, we show that the adaptation to high-temperature and acidic environments of these lineages was driven by different strategies, with the prior relying more on innate genomic characteristics and the latter more on the acquisition of genes associated with acid tolerance. Finally, by reconstruction of ancestral states of OGT and pI, we showed the potential evolutionary process of Parvarchaeales-related lineages with regard to the shift from a high-temperature environment of their common ancestors to low-temperature (potentially acidic) environments.

COVID-19 vaccination resistance has become a major challenge to prevent global SARS-CoV-2 transmission. Here we report that the vaccination coverage rate is inversely correlated to the mutation frequency of the full genome (R2=0.878) and spike gene (R2=0.829) of SARS-CoV-2 delta variants in 16 countries, suggesting that full vaccination against COVID-19, with other mitigation strategies, is critical to suppress emergent mutations. Neutrality analysis of DH and Zengs E tests suggested that directional selection was the major driving force of delta variant evolution. To eliminate the homogenous effects (population expansion, selective sweep etc.), the synonymous (Dsyn) and nonsynonymous (Dnonsyn) polymorphisms of the delta variant spike gene were estimated with Tajimas D statistic. Both D ratio (Dnonsyn/Dsyn) and {Delta}D (Dsyn-Dnonsyn) have positive correlation with the full vaccination rate (R2= 0.723 and 0.505, respectively) in 19 countries, indicating that purifying selection pressure of SARS-CoV-2 spike gene increased as the vaccination coverage rate increased. Taken together, our data suggests that vaccination plays an important role in the purifying selection force of spike protein of SARS-CoV-2 delta variants.

The outbreak of coronavirus disease 2019 (COVID-19) which originated in Wuhan, China, constitutes a public health emergency of international concern with a very high risk of spread and impact at the global level. We developed data-driven susceptible-exposed-infectious-quarantine-recovered (SEIQR) models to simulate the epidemic with the interventions of social distancing and epicenter lockdown. Population migration data combined with officially reported data were used to estimate model parameters, and then calculated the daily exported infected individuals by estimating the daily infected ratio and daily susceptible population size. As of Jan 01, 2020, the estimated initial number of latently infected individuals was 380.1 (95%-CI: 379.8[~]381.0). With 30 days of substantial social distancing, the reproductive number in Wuhan and Hubei was reduced from 2.2 (95%-CI: 1.4[~]3.9) to 1.58 (95%-CI: 1.34[~]2.07), and in other provinces from 2.56 (95%-CI: 2.43[~]2.63) to 1.65 (95%-CI: 1.56[~]1.76). We found that earlier intervention of social distancing could significantly limit the epidemic in mainland China. The number of infections could be reduced up to 98.9%, and the number of deaths could be reduced by up to 99.3% as of Feb 23, 2020. However, earlier epicenter lockdown would partially neutralize this favorable effect. Because it would cause in situ deteriorating, which overwhelms the improvement out of the epicenter. To minimize the epidemic size and death, stepwise implementation of social distancing in the epicenter city first, then in the province, and later the whole nation without the epicenter lockdown would be practical and cost-effective.

CRISPR/Cas9 mediated precise gene editing requires homology-directed repair (HDR), which occurs less frequently than non-homologous end-joining (NHEJ) including the canonical NHEJ and alternative NHEJ (Alt-EJ) in mammalian cells, especially in CHO cells that inherent resist HDR. To solve the above hurdle, here we for the first time show that knockout the DNA polymerase {theta} (POL{theta}), which is essential for Alt-EJ, significantly increases the knock-in efficiency by nearly forty-fold in CHO cells via eGFP reporter system and does not affect the normal growth and proliferation of cells. Meanwhile, even when transfecting simple circular, without negative element homologous template DNA donor and CRISPR/Cas9 plasmid to two different genomic sites, the knock-in rate of 4kb donor integration can still reach a mean of over 80% (29/36) and 2.7% (1/36) of the selected cell colonies in POLQ-/- CHO cells, however, no positive knock-in cell colonies was obtained in wild-type CHO cells which respectively selected 62 cell colonies and 36 cell colonies. Furthermore, we show that POLQ promotes random integration in CHO cells. Finally, RNA-sequence analysis reveals not significant altered DNA repair, metabolism, apoptosis, and cell cycle in POLQ-/- cells. These findings open a new target gene POLQ to overcome bottlenecks of the precision genome editing.

Amphioxus is considered the best-known living proxy to the chordate ancestor and an irreplaceable model organism for evolutionary studies of chordates and deuterostomes. In this study, a high-quality genome of the Beihai amphioxus, Branchiostoma belcheri beihai, was de novo assembled and annotated. Within four amphioxus genomes, twenty-eight groups of gene novelties were identified, revealing new genes that lack homologs in non-deuterostome metazoa, but share unexpectedly high similarities with those from non-metazoan species. These gene innovation events have played roles in amphioxus adaptations, including innate immunity responses, glycolysis, and regulation of calcium balance. The gene novelties related to innate immunity, such as a group of lipoxygenases and a DEAD-box helicase, boosted amphioxus immune responses. The novel genes for alcohol dehydrogenase and ferredoxin could aid in the glycolysis of amphioxus. A proximally arrayed cluster of EF-hand calcium-binding protein genes were identified to resemble those of bacteria. The copy number of this gene cluster was negatively correlated to the sea salinity of the collection region, suggesting that it may enhance their survival at different calcium concentrations. This comprehensive study collectively reveals insights into adaptive evolution of cephalochordates and provides valuable resources for research on early evolution of deuterostomes.

Acinetobacter baumannii is an emergency pathogenic bacterium for its multidrug-resistance and high mortality rates after infection. In-depth genetic analysis of A. baumannii virulence and drug-resistant genes is highly desirable. Existing methods for genetic manipulation of A. baumannii mainly rely on the use of antibiotic as the selectable marker, and the sacB/sucrose as the counter-selectable marker, which is inconvenient and inappropriate for all research of A. baumannii. Based on the highly conserved pyrF gene and its conserved 500bp-flanking sequence, we quickly and easily generated the pyrF-deleted mutants as the uracil auxotrophic host strain in three model strains and 11 clinical strains. The pyrF-carried vectors constructed for gene editing with pyrF/5-FOA as the counter-selection were conveniently and time-saving in these pyrF-deleted mutants. Utilizing the pyrF-based genetic manipulation system, we easily and efficiently modified the cas gene and CRISPR sequence of I-F CRISPR-Cas system in A. baumannii AYE, and detected the CRISPR interference and adaptation in these mutants. In summary, the pyrF-based genetic manipulation system could be broadly applicable used for efficiently maker-less gene editing in most A. baumannii strains.

The outbreak of the novel coronavirus disease, COVID-19, originating from Wuhan, China in early December, has infected more than 70,000 people in China and other countries and has caused more than 2,000 deaths. As the disease continues to spread, the biomedical society urgently began identifying effective approaches to prevent further outbreaks. Through rigorous epidemiological analysis, we characterized the fast transmission of COVID-19 with a basic reproductive number 5.6 and proved a sole zoonotic source to originate in Wuhan. No changes in transmission have been noted across generations. By evaluating different control strategies through predictive modeling and Monte carlo simulations, a comprehensive quarantine in hospitals and quarantine stations has been found to be the most effective approach. Government action to immediately enforce this quarantine is highly recommended.

It has been challenging to characterize the lineage relationships among cells in vertebrates, which comprise a great number of cells. Fortunately, recent progress has been made by combining the CRISPR barcoding system with single-cell sequencing technologies to provide an unprecedented opportunity to track lineage at single-cell resolution. However, due to errors and/or dropouts introduced by amplification and sequencing, reconstruction of accurate lineage relationships in complex organisms remains a challenge. Thus, improvements in both experimental design and computational analysis are necessary for lineage inference. In this study, we employed single-cell Lineage tracing On Endogenous Scarring Sites (scLOESS), a lineage recording strategy based on the CRISPR-Cas9 system, to trace cell fate commitments for zebrafish larvae. With rigorous quality control, we demonstrated that lineage commitments of complex organisms could be inferred from a limited number of barcoding sites. Together with cell-type characterization, our method could homogenously recover lineage information. In combination with the cell-type and lineage information, we depicted the development histories for germ layers as well as cell types. Furthermore, when combined with trajectory analysis, our methods could capture and resolve the ongoing lineage commitment events to gain further biological insights into later development and differentiation in complex organisms.

Targeted insertion of exogenous sequences to genomes is useful for therapeutics and biological research. While CRISPR/Cas technologies have been very efficient in gene knockouts by double-strand breaks (DSBs) followed by indel formation through non-homologous end-joining (NHEJ) repair pathway, the precise introduction of new sequences mainly rely on inefficient homology directed repair (HDR) pathways following Cas9-induced DSBs and are restricted to dividing cells. The recent invention of Prime Editing allows short sequences to be precisely inserted at target sites without DSBs. Here, we combine Prime Editing and sequence-specific recombinases and integrases to insert kilobase sequences directionally at target sites. This technique, called PRIMAS for Prime editing, Recombinase, Integrase-mediated Addition of Sequence, will expand our genome editing toolbox for targeted insertion of long sequences up to kilobases and beyond.

The CRISPR/Cas9 platform holds promise for modifying fish traits of interest as a precise and versatile tool for genome manipulation. To reduce introgression of transgene and control reproduction, catfish species have been studied for upscaled disease resistance and intervening of reproduction to lower the potential environmental risks of introgression of escapees as transgenic animals. Taking advantage of the CRISPR/Cas9-mediated system, we succeeded in integrating the cathelicidin gene from an alligator (Alligator sinensis; As-Cath) into the target luteinizing hormone (LH) locus of channel catfish (Ictalurus punctatus) using two delivery systems assisted by double-stranded DNA (dsDNA) and single-stranded oligodeoxynucleotides (ssODNs), respectively. In this study, high knock-in (KI) efficiency (22.38%, 64/286) but low on-target was achieved using the ssODN strategy, whereas adopting a dsDNA as the donor template led to an efficient on-target KI (10.80%, 23/213). On-target KI of As-Cath was instrumental in establishing the LH knockout (LH-_As-Cath+) catfish line, which displayed heightened disease resistance and reduced fecundity compared to the wild-type sibling fish. Furthermore, implanting with HCG and LHRHa can restore the fecundity, spawnability and hatchability of the new transgenic fish line. Overall, we replaced the LH gene with an alligator cathelicidin transgene and then administered hormone therapy to gain complete reproductive control of disease-resistant transgenic catfish in an environmentally sound manner. This strategy not only effectively improves the consumer-valued traits, but also guards against genetic contamination. This is a breakthrough in aquaculture genetics to confine fish reproduction and prevent the establishment of transgenic or domestic genotypes in the natural environment.

Humanized rodent models, especially humanization of genetic/genomic components involved in immunity have significantly advanced our understanding of human immune system. Here, we utilized trans-chromosomic (Tc) technology to generate a TC-mAb rat model that stably harbors a mouse artificial chromosome carrying full-length human immunoglobulin (Ig) heavy and kappa light chain genes (IGHK-NAC) in a rat Ig knockout background. In contrast with TC-mAb mice, serum human IgG concentration was found higher than IgM. Number of lymphocytes was recovered, and B cell population in the spleen was normal. Remarkably, repertoire analysis revealed similarities between the model and human PBMCs; somatic hypermutation and class-switch recombination also more closely resembled humans. Furthermore, immunization resulted in generation of antigen-specific human antibodies. Collectively, our strategy to generate both rat and mouse models through introduction of the identical IGHK-NAC offers unprecedented opportunities to comprehensively evaluate genomic regulation and its outcomes associated with genomic sequences and host-derived protein factors.